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Abstract In this paper we are interested in Euler-type sums with products of harmonic num- 
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given. The given representations are new. 
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1 Introduction and Preliminaries 


1.1 Harmonic numbers and Euler sums 


The generalized nth-harmonic numbers of order p are given by ( [23]) 


1 
HP :=Y 3, pen, (1.1) 
j=1 


where the H,, = HP is the classical harmonic number ( [1]) and the empty sum He ) is con- 
ventionally understood to be zero. 
In response to a letter from Goldbach in 1742, Euler considered sums of the form (see 


Berndt [4]) 
ES aw) 
0 
n=1 


where p, q are positive integers with q > 2, and w := p+q denotes the weight of sums Sy y. These 
sums are called the linear Euler sums today. Euler discovered that in the cases p = 1, p = q 
and p + q is less than 7 or when p + q is odd and less than 13, the linear sums have evaluations 
in terms of zeta values, i.e., the values of the Riemann zeta function ( [2]) 


at the positive integer arguments. Moreover, he conjectured that the double linear sums would 
be reducible to zeta values when p+q is odd, and even gave what he hoped to obtain the general 
formula. In 1995, Borweins and Girgensohn [12] proved conjecture and formula, and in 1994, 
Bailey, Borwein and Girgensohn [3] conjectured that the linear sums S,, y, when p+q > 7 and 
p+ q is even, are not reducible. Hence, the linear sums S,,y can be evaluated in terms of zeta 
values in the following cases: p = 1,p = q,p + q odd and p+q= 4,6 with q > 2. 

Similarly, the nonlinear Euler sums are the infinite sums whose general term is a product 
of harmonic numbers of index n and a power of n”?. Let rm := (71,...,Ty) be a partition of 
integer p and p=71 +-+::** +7 with 11 < 79 <-:+: < Tp. The classical nonlinear Euler sum of 
index rr, q is defined as follows (see [23]) 


00 


AO gl) 0% ar 


m3 , 922 (1.2) 


Srir2:T1:9 > 
n=1 


where the quantity w := 71 +-::* +7 + is called the weight and the quantity k is called the 
degree. As usual, repeated summands in partitions are indicated by powers, so that for instance 


00 2) (4) 
EH 
S12234,q = 91122249 = > a 
n=1 


The nonlinear Euler sums, i.e., Sr y with 7 having two or more parts, are more complicated. 
Such sums were already considered in [3,23, 43,46, 47, 49,51,52]. In [23], Flajolet and Salvy gave 
an algorithm for reducing Sy; a, to linear Euler sums when 71 + 712 + q is even and 71,72,q > 1 


(see Theorem 4.2 in the reference [23]. In [46,47], the first author jointly with Li, Yan and Shi 
proved that all quadratic Euler sums of the form 


00 Pp 
H,H, 
Sita = E á—Á (p 922, p+9<10) 


can be reduced to polynomials in zeta values and linear sums. Moreover, in the recent papers 
(43, 49], the first author [49] shown that all Euler sums of weight < 8 are reducible to Q-linear 
combinations of single zeta monomials with the addition of (S26) for weight 8. For weight 9, 
all Euler sums of the form Sr;...rr,, With q € (4,5,6,7) are expressible polynomially in terms 
of zeta values. Wang and Lyu [43] shown that all Euler sums of weight eight are reducible to 
linear sums, and proved that all Euler sums of weight nine are reducible to zeta values. In the 
most recent papers [44,51,52], the first author has carried out fruitful cooperation with Wang 
proved that all Euler sums of weight ten can be expressed as a rational linear combination of 
¿010,265.80 0865), 4(3C(4),0(2) 26 and S2g. All Euler sums of weight 
11 are reducible to 926, 618,2 and zeta values. 

Similarly, in [23], Flajolet and Salvy introduced the following alternating harmonic numbers 


n aqi l a ml 

HP =yY Ea and H,:=H(P, pen. (1.3) 
: j 
j=1 


Clearly, we have the following relations for harmonic number and alternating harmonic number 
q) — yq) A ( (») (») ) 
1 Paya 7 29 Hija + H—/2) 222 


He = Hn — Hin/2)> 


where [x] denotes the greatest integer less than or equal to z. Ñ 
The alternating harmonic number converges to the alternating Riemann zeta value C(k): 


lim AY =((k), keN, 


n—>00 


where the alternating Riemann zeta function is defined by (for more details, see [23]) 


o (_qynm-1 
a y SA Rs) > 1: 
n=1 


ns ” 


The generalized alternating Euler sums involving harmonic numbers and alternating harmonic 
numbers are defined by 


oo prlpr) (m1) q(ri) (rr) 
ba a AA (14) 
O» ra 1 de 
i=1  j=1 
o ae) em aer) en ar) me 
O) a e E A 15) 
T » Il 75.4 E n 


where the (py,...,p1) € (NY, (r1,...,7p) € (NY with p; < --- < p¡ and r¡ < --- < rp and the 
quantity w:=p1+::**+p +11 +-**+Tr£+q being called the weight and the quantity l+k being 


the degree, l > 0 and k > 1. Since repeated summands in partitions are indicated by powers, 
we denote, for example, the sums 


a SHA A Ama aa 
132123,5 — yd ——= 7 3 123133,4 — A CN: 
n=1 n=1 


It has been discovered in the course of the years that many alternating Euler sums admit 
expressions involving finitely the zeta values, ln2 and polylogarithms. For example, in [23], 
Flajolet and Salvy gave explicit reductions to zeta values and logarithm for all linear sums 


co aw 00 q») 00 aw de 00 aw Pm 
Spy =D BRE y ao Pas y na 1)" and Sp q = y aL) 
n=1 n=1 n=1 n=1 


when p + q is an odd weight. The evaluation of the above four linear sums in terms of values 
of the Riemann zeta function and polylogarithm function at positive integers is known when 
(p, q) = (1,3), (2,2), or p+q is odd (see [12,20,23,45]). Here are a few explicit evaluations 


7 5 

S13= ze (3) In 2 — 104), 

Ei = ADE a + Lc (9)Im22 2 pto — a (3) In 2 
1,3 — 14 A 9 12 Fl , 

S13= sE (4) + E (2) In2 — In = Mi (12%, 

S23 = 604) + 4Li (1/2) + 1c(3) In 2 — ((2)In?(2) + Ina 


Here Li,(x) denotes the polylogarithm function, which is defined for |+| < 1 by 


n=1 ds 
Furthermore, we consider the nested sum 
0,8) n n k 
y z 
y — |, 17,y€el|-1,1), mpenN. 
m P > y) y) ? e) 
n=1 ds (2 h 


By taking the sum over complementary pairs of summation indices, we obtain a simple reflection 


formula 
oo n 


DD) Fi) Pro tin to + Len. 
=1 k=1 n=1 k=1 


n 


Hence, we readily deduce the following four results of linear sums 


13 19 3 

535 == (4), S3i= 10D = ze (3) In 2, 

Saa = E (1) AL (1/2) LO) né2— ina Pa le (3) In 2, 
(1 3 1 Ñ 1 

S3 1 =2Lia (3) E (3) In 2 + pre? e (4) — as (2) In?2. 


However, in many other cases we are not able to obtain a formula for the Euler sum constant 
explicitly in terms of values of the Riemann zeta, logarithm and polylogarithm functions, but 
we are able to obtain relations involving two of more Euler sum constants of the same weight. 
In [48], the first author jointly with Yang and Zhang proved that the alternating quadratic Euler 
sums 


[0.0] [0.e] => [e.e] => 
He Ae HB? Y An Hn =1 

Y ED Sra 2120") $135 = » a E 

n=1 n=1 n=1 

[0.0] => [e.el 00 

H, H. HB? HB? 

Sii om = o —=2 and A 

2 n2m 'm 2 n2m a (2n al 1) m 


are reducible to polynomials in zeta values and to linear sums, and are able to give explicit 
values for certain of these sums in terms of the Riemann zeta function and the polylogarithm 
function, here m is a positive integer. In [3], Bailey, Borwein and Girgensohn considered sums 


Y Hin+ 1)”, Y Hln+ 1) PD" and Y HU DEDO, (p+as5) 
n=1 


n=1 


and obtained a number of experimental identities using the experimental method, where p > 2 
and q > 1 are positive integers. 

Euler sums may be studied through a profusion of methods: combinatorial, analytic and 
al- gebraic. There are many other researches on Euler sums and Euler type sums. Some related 
results for Euler sums may be seen in the works of [13, 14,16, 17, 24, 32,33, 39-41, 50,53]. 


1.2 Multiple harmonic (star) numbers and multiple zeta (star) values 


For m € N, $ := (s1,82,...,8m) € (N)”, and a non-negative integer n, the multiple harmonic 
number (MHN for short) and the multiple harmonic star number (MHSN for short) are defined 
by, respectively 


1 
Ch (its) ¡== ME (1.6) 
Lc km 
1<km<:-<kj] <n 
1 
Gh (81,82, --- 18m) := TS (1.7) 


1<km<--<k¡<n 


The integers m and w := |S| := es s¡ are called the depth and the weight of a multiple 
1= 


harmonic number or multiple harmonic star number. By convention, we put (. (S) = 0 if 
n < m, and (n (0) = (0) = 1. By [s1,8s2,...,sj¡)”. we denote the sequence of depth mj with 
m repetitions of (s1,52,...,S7). 


If s; > 1, then the limits as n > oo exist, and are the well-known multiple zeta and multiple 
star-zeta values ((51,892,...,8m) and C* (51,89,...,Sm) respectively. We put a bar on top of 
si (j = 1,:-:k) if there is a sign (—1)”% appearing in the denominator on the right. For 
example, 


(=p A+ 
Gn (31, 82, 83, --+>8m—1,8m) = E ET (1.8) 
1 m 


1<km<--<k1<n 


k2+Km 

OS - (1) 
€ (81432, 838-143) = 8... fém' 
1<km<-"<k¡<n 1 ls 


(1.9) 


The sums of types (1.8) and (1.9) are called the alternating harmonic numbers and alternating 
harmonic star numbers, respectively. For alternating MHNs or MHSNs, we have the similar 
limit cases. For example 


Ep 
C(S5L3% <<: Em) i= lim Ón (31,82,+=<.8m) = 2 RE Em? (1.10) 
1<km<"<k] 1 mu 
¡Ent 
*k = = i= . k — = a 
EA Lim Cm (81,580: .8m) = y TOS (1.11) 


1<km<-"<kis<n 1 


The sums of types (1.10) and (1.11) (one of more the s; barred) are called the alternating 
multiple zeta values and alternating multiple zeta star values, respectively. 
If the first index is barred, then the limits as n > oo exist for all positive-integer indices, e.g., 


¿(1 = lim Cn (1) =-In2. 


Multiple zeta values are certainly interesting and important. The number C(6,2) appeared in 
the quantum field theory literature in 1986 [15], well before the the phrase multiple zeta values 
had been coined (Multiple zeta values were introduced by the Hoffman [26] and Zagier [55] in 
the 1990s). The alternating multiple zeta values have been intensely studied since the early 
1990s, and extensive tables of them were compiled, first by the group at Lille [6], and later as 
part of the Multiple Zeta Value Data Mine [7]. 

In the past several bases have been considered for both the multiple zeta values and al- 
ternating multiple zeta values. The vector space of multiple zeta values can be constructed 
allowing basis elements, which contain besides the zeta values the index of which is a Lyndon 
word products of this type of zeta values of lower weight. One basis of this kind is (see [7,34,36]) 


w=2, (2), 

w=3, 6), 

w=4, C(4), 

w=5, (5), E(2)08), 

w=6, C(6), (3), 

w=7, C(7), C(2)05), (8B)C(4), 

w=8, C(8), E(3)C(5), C(2)C(3), (65,3), 

w=9, C(9), E(2)E(7), E(3)C(6), C()C(5), (3), 


w=10, (10), E(3)C(7), (5), C(2)6(3)C(5), E(4)%(3), C(2)C(5, 3), C(T, 3), etc. 


Rational relations among multiple zeta values are known through weight (sum of the indices) 
22 and tabulated in the Multiple Zeta Value Data Mine (henceforth MZVDM) [7]. 

A basis for the Q-vector space spanned by the set of alternating multiple zeta values for 1 < w< 6 
should be given 


) 
). (2) 2, In*2, 

E (1/2), C(4), C(3) In 2, c(2) In? 2, In? 2, 

w=5, Lis(1/2), C(5), Li4(1/2)In2, C(2)C(3), C(4)In2, C(3)In*2, C(2)In*2, In?2, 

= 1 ( ), (6), C(5, D) Lis(1/2) In 2, (5) In 2, Li4(1/2)C(2), Li1(1/2) In? 2, (3), 
1:(3)In2, C(4)In?2, ((3)In92, ((2)In*2, In* 2, etc. 

Rational relations among alternating multiple zeta values are also tabulated in the MZVDM 
(through weight 12 [7)). 

The relations between multiple zeta (star) values and the values of the Riemann zeta function 


and polylogarithm function have been studied by many authors. For explicit evaluations of the 
multiple zeta (star) values or related alternating sums, readers may consult [9-11,22,27,29,32,35, 


42,49,55-57]. For example, Zagier [56] proved that the multiple zeta star values (* (ey, 701 2) 
and multiple zeta values ( (23, a, er), where a, b € No := NU[O] = [0,1,2,...), are reducible 


to polynomials in zeta values, and gave explicit formulae. Hessami Pilehrood et al. [36,37] and 
Li [30] provide three new proofs of Zagier”s formula for (* (ey, E (y) based on a finite identity 
for partial sums of the zeta-star series and hypergeometric series computations, respectively. 


1.3 Relations between Euler sums and multiple zeta values 


From the definitions of harmonic numbers and multiple harmonic numbers, we can find that 


=((p) and HP =-((p), pen. 


According to the rules of the “harmonic algebra” or “stufle product”, it is obvious that the 
products of any number of harmonic numbers and alternating harmonic number can be expressed 
in terms of (alternating) multiple harmonic numbers. For example, 


q e ¿0)+24 11,1), 
HH? = Gn (1) Gn (3) = Gn (1,3) + Gn (3, 1) + Gn (4), 
HHa == 1 ) ón Me (10=(1,0=4 (1,1): 


Note that when different indices are combined, the combination of a barred and unbarred index 
is barred, while the combination of two barred indices is unbarred. 
The series 


2 , 7 Se. ¿7 . 
y nl . E). ij EN 
can be written as the sum of multiple zeta value 
Clp,t1, 12, > 16) HC l(p+ 1,42): ,Íp)- 


Series of the form 


o Cn (11,12, *** ,1y) Z Cn (01, 19,** 06) n 
a E E EA E ES E Y 
n? nP 

n=1 n=1 


with one of more of the 1; barred, can be written as a sum of alternating multiple zeta values. 
It is clear that every (alternating) Euler sum of weight w and degree n is clearly a Q-linear 
combination of (alternating) multiple zeta values or (alternating) multiple zeta star values of 
weight w and depth less than or equal to n+1. In other words, (alternating) multiple zeta (star) 
values are “atomic” quantities into which (alternating) Euler sums decompose. For example, 
St 2 = 6 (2, 2) =E 26 (2, L, il = 4 (4) == 26 (3, 1) > 

51314 = € (4,1,3) + € (4,4) + € (4,3,1) + € (5,3) + € (8) + € (7,1), 

Sis ==((,2) -£(2,1,1) -((2,1,1) -( (4) = C (3, ip c (3,1). 
Hence, based on the above, we obtain the following two conclusion: 

1. All Euler sums of weight € [3,4,5,6,7,9) can be expressed as a rational linear combi- 
nation of products of zeta values. For weight 8, all such sums are the sum of a polynomial in 
zeta values and a rational multiple of C(5,3). For weight 10, all such sums are reducible to zeta 
values, ((5,3) and C(7,3). (All explicit formulas can be found in [23, 43, 44, 46, 47, 49)) 

2. All alternating Euler sums of weight < 5 can be expressed in terms of zeta values, 
polylogarithms and In 2. (Apart of explicit formulas can be found in (3, 12,33, 45, 48]) 

In this paper, we will give almost all explicit evaluations of alternating Euler sums with 
weight less or equal to 5. 

The evaluation of Euler sums and multiple zeta values has been useful in various areas 
of theoretical physics, including in support of Feynman diagram calculations and in resolving 
open questions on Feynman diagram contributions and relations among special functions [18], 
including the dilogarithm, Clausen function, and generalized hypergeometric function [19]. An 
array Of Euler sums and multiple zeta values is required in calculations of high energy physics. 
These quantities appear for instance in developing the scattering theory of massless quantum 


electrodynamics [8]. Harmonic number sums also occur in computer science in the efficiency 
analysis of algorithms [3]. 


1.4 Purpose and plan 


The purpose of this paper is to evaluate (alternating) Euler sums and related sums involving 
(alternating) harmonic numbers. 

The remainder of this paper is organized as follows. 

In the second section we evaluate some Euler related sums which involve harmonic num- 
bers (or alternating harmonic numbers), Stirling numbers and Bell numbers. Then we use the 
evaluations obtained to establish many identities involving two or more Euler sums of the same 
weight, which can be expressed in terms of Riemann zeta values and polylogarithms. 

In the third section we use certain integral representations to evaluate several series with 
harmonic numbers and alternating harmonic numbers. Moreover, we establish many relations 
involving two or more Euler sums of the same weight. Using these relations obtained we give 
many closed form of alternating Euler sums with weight < 5. Furthermore, we also consider the 
following type of Euler-type sums 


-Y Cul MTL) “Ca (Mi Tp) 


A k N 
(n +1) MEE TIT”, P)K,TEN, 


Sm1ma...mplF, k; Ll). Up, L 


and provide some explicit evaluation in a closed form in terms of zeta values and harmonic 
numbers. 


2 Series involving Stirling numbers, harmonic numbers and Bell 
numbers 


In this section, we will establish some explicit relationships that involve Euler sums and zeta 
values. 


2.1 The Bell polynomials, Stirling numbers and Bell numbers 
The exponential complete Bell polynomials Y,, are defined by (see section 3.3 in the [21] and 
section 2.8 in the [38]). 


qe E 
exp Nam =1+) Y, (51, La, => Bj) Pr (2.1) 


ml 
m>1 k>1 


Then Yo (-) = 1 and 
k! 21 Y 41 (Lay 22 LpN 5h 
Amaia] Y a A 
c1+2c2+-"+kcp=k 


Additionally, by Eq. (2.44) of [38], the polynomials Yz.(-) satisfy the recurrence relation 


k-1 
k-—1 
Ye (aia 2) Y ( E ) ed (aaa 23) ken. (2.3) 
j=0 


The (unsigned) Stirling number of the first kind s(n, k) is defined by [21] 
| a pe 
nix (1+x2x) (1+ 5) (1+2) pd ; (2.4) 


with s(n, k) =0,ifn < kand s(n,0) = s(0,n) = 0, s(0,0) = 1, or equivalently, by the generating 
function: 


iO pS sei (2.5) 


n=k 


Moreover, the (unsigned) Stirling numbers s(n, k) of the first kind satisfy a recurrence relation 
in the form 


s(n,k)=s(n-=1,k-1)+(n-1)s(n—1,k). (2.6) 


By the definition of s(n, k), we see that we may rewrite (2.4) as 


n 


a z 
s(n+1,k+1)]a* =n!exp in (142) 


k=0 j=1 


n 00 k 
= nlexp A 


j=1 k=1 


00 (k) 
ono DE) a (2.7) 


k=1 
Letting 2 = (-1)%7 (% — Das in (2.1) and combining (2.7), the following relation holds: 
_ mM re) k-1 rre) 
s(n+1,k+1) = mel Ha, UE) (EDO Ak DIE). (2.8) 


Therefore, we know that s(n, k) is a rational linear combination of products of harmonic numbers. 
Then by (2.3), we have the recurrence 


k-1 
s(n+1,k+1)= pl DEFIgeDs(n+1,j+1), ken, (2.9) 
22 
which further yields 
s(n+1,1) =n!, s(n+1,2) =n!H,, s(n+1,3) = ” [m3 - Hp | A 


I 
s(n+1,4) = 5 Es 34. HD 4280, 


s(n+1,5) = "mi CHO — 62H +3(40 y +8H, H6 1. 


24 
Throughout this article we define a Bell number (Y¿(n)) by the Bell polynomials 
Y, (n) = Y, (4. LH... (k- 1110) (2.10) 
Thus, according to the recurrence (2.3), we obtain 
BL 
Yom) =1, Y) => (577) DD (a), 
j=0 


from which, we know that Y; (n) is a rational linear combination of products of harmonic 
numbers, and we have 


Y (nm) = E (nm) = Hi + HO, (1) = Ei +SE HPF 2ED, 
Y, (n) = Hi +8H,H% + 68240 + 3407? +68, 


2.2 Some lemmas 


The following lemmas will be useful in the development of the main theorems. 


Lemma 2.1 (/49)) For integers m > 0 and n > 1, then the following integral identity holds: 


Ie Mat "d= 27 5 a O mayra, zx € (0,1). (211) 


0 


10 


Lemma 2.2 For positive integers m and nonnegative integers k, 


1 
ln” (1 — 2) ln* 
W (m, k) = PA (2.12) 
z 
0 
then 
W (m, k) = (-1)"+miktc (1 +2, pq (2.13) 
Proof. Expanding out ln”"(1 — x) and using the lemma 2.1, 
1 
de 00 00 00 dd Hr 
mr) NY as 
i1=142=1 Jas tm 


q DR! 


até nc 01120 imlt1 ++ +ém 
By Lemma 4.3 in the [26], 


o PP 


ii=1i9=1  im=1?1?2* Ep (1 + im) n1>n2>-->Nm >1 


and the conclusion follows. 
In [9], Borwein, Bradley and Broadhurst gave the generating function 


> C (m + a E ma"? 1 —exp (Sc | (2.14) 


n,m=1 


which implies that for any m, n € N, the multiple zeta value € (m +1, ay) can be repre- 


sented as a polynomial of zeta values with rational coefficients, and we have the duality formula 
E (n +1, na =C (m did, y) 
In particular, one can find explicit formulas for small weights. 
(2, (17%) =C(m+2), 
c(8,(1)») = mz Ea)= E cn 


In [46], the first author jointly with Yan and Shi gave the following recurrence relation of W (m, k) 


m-—1k-1 
W/O) (G) ran), (15) 


i=1 j=0 
where 4 (z) denotes the digamma function (or called Psi function ) which is defined as the 
logarithmic derivative of the well known gamma function: 
Pda 
(2) 
From (2.14) or (2.15), we see that the values of integrals W(m, k) can be expressed as a rational 
linear combination of products of zeta values. 


— 


e 


Y (2) :=, (nT (2) = 


11 


Lemma 2.3 For positive integers m and nonnegative integers k, 


aja TS (1+2x) In" EN (2.16) 
0 
then 
V (m,k) = (DIH mk (+2, py) (2.17) 


Proof. The proof of Lemma 2.3 is similar to the proof of Lemma 2.2. 


Lemma 2.4 Ifm> 1 is a integer and z € [0,1] , then we have 


pa, a) 


z m-+1 1+2 
0 
sí l 
1 m=-j+1 1 1 
y) a == (2.18) 
ce (m-—j+1)! l+2 


Proof. We note that the integral in (2.18), which can be rewritten as 


z : de (1+2)71 1 
m 1 = a mE ap =1 Pm m 
pe o ] — a A ] du 
dl t-1 u—u? 
0 1 dl 
(+2) (1+2) e 
In"u aa 
1 m+1 
(1) 1 dos 1] 4 
1 1 
(1+2) e 
1 ay 
=—— A 1 yr J 2.19 
re (1 2) + (1) da (249) 
1 
Then with help of formula (2.11) we may deduce the desired result. 
Lemma 2.5 (see [46/) For integers n > 1 and k > 0, then 
/ Y; 
Jenna at = (ap (2.20) 
n 
0 


where the Bell number Yi (n) is defined by (2.10), and Yy (n) = k!ICf (694) ; 


Lemma 2.6 (/5,25/) For positive integer n, m and any suitable sequences a; EC (¡=1,2...,n), 


lo =D) Amín)t”, Ao(n)=1, 


i=1 m=0 


][ ( + at) = 7 Am(n)t”, Ap(n)=1. 


i=1 m=0 
then 
Asin)j:= Ak *** Oli (2.21) 
1<km<-<k¡ <n 
Ar (n) := y lez *** O (2.22) 
1<km<-"<k1<n 
] -1 =1 
Here |t| < min ( ja] ¡e al ) : 


2.3 Main results and proofs 


Theorem 2.7 For positive integer m, then the following symmetric-function identities hold: 


1 m-—1 Epa m-—1 
him == hiPm=i and em = ———— (1) eiPm-i, (2.23) 
m m : 
i=0 i=0 
where hm = hm[(21,22,*** ,Tn) and €m = €m[11,T2,*** ,Tn) are respectively the complete and 


elementary symmetric functions of degree m, and Pm = 1" +13 +---+a7 is the mth power 
sum (for details introductions, see [28/, the notation follows [31/). 
Proof. From Hoffman's paper [28], we know that 


[0] n 


Be =TH (1+tx;) and HA(t => yl a |] —- 
i=1 


j=0 i=1 


are the respective generating functions of the elementary and complete symmetric functions. 
Note that we can rewrite E(t) and H(t) as 


EA OL 
E (t) = exp >» — ot and H(t) =exp EN pont ; 
k=1 k=1 
Here |t| < min (Ja ]7l, +. O Letting 2, = (-1)**(k — 1)!pz and (k — 1)!p; in (2.1), 


then comparing the coefficients of t*, we can find that 


Y, (os, —1lpa, >>> (1) (k — 1)'px) 
€ = and h;= 


Y, (pr, 1!p2, ER (k —> 1)!px) 
kl : 


k! 


Using the recurrence relation (1.3), by a simple calculation, we obtain the desired results. 
(Theorem 2.7 can also be found in Entry 12(a) of Chapter 14 in Ramanujan's second notebook 
[5]). Hence, the elementary symmetric function ez and the complete symmetric function hz can 
be expressed in terms of the power sumas, i.e., (see [28]) 


Ek = Ek (p1,pa,-*: , Pr) and hr ES Q (p1, pa," , Pr) - 
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Moreover, from the generating function E(t) and H(t), we have 


[e.e) n [0.0] 
1 
O Em (21,*** ,En)t” = [10 ++.) = ¡EN y Em (L1,*** Tn, Tn 1) 8" 
m=0 4=1 m=U 
[oe m 
3 (Eco das 0) en 
m=0 Mk=0 
Y ee 1 1 a 
y Pon (Di Tp, Tp +1) 1" > 1L71=2x, = Tfn y domi (ij Dr 
im=0 i=1 m=0 
[0 m 
E Estatutos m]e0 
m=0 Xk=0 


Then by considering the coefficient of t””, we can get the following relations 


m 
Em (“a 2 ¡Ea) a y (DE em (Dis ES Brida 
k=0 
m 
an (La, poÑos ni) = ala ESP GE El 
k=0 


The polynomials P,, and Q,, can be written as determinants (see [29]) 


y y 2 +. 0 
AAA A O: pos : 
Vit UD 31 1 


Yn Yn-1Yn-2***  Y1 
and 
Yi —1 0 0 
Ya Yi —2 0 
Yn-—1 Yn-2 Yn-3 ::: =N =+ 1 
Yn Yn-1 Yn-2*** Y1 


Here are a few explicit evaluations; some of them were given previously by Hoffman and Ohno. 


Y — Y2 
Pi (y) =y1, Poly, ya) = a" 


3 
yí — 3y1Y2 + 2y3 
Ps (y, ya, ya) = ARA 6 ; 


4 2 2 
yí — 6yíy2 + 8Y1Y3 + 3y2 — Oya 
Pa (yr, ya, ya, ya) = =— AAA E 94 = ; 


yí — 10y¿yo + 20yÍy3 + 15y143 — 30y1 ya — 20y2y3 + 24ys 
Ps (Y1, Y2) Y3, Ya, Y5) = a --- 
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and 


yi +Y 
Qi(u)=v4, Qoíyi ya) = 2, 


2 
y] +3Yy1iYo + 2y3 
Q3 (yr, Y2, Yy3) = CI 
yi + 6yiyo + 8y1Y3 + 3y3 + 6ya 
a (Yi, Y2, Ya, Ya) = ERE SRT: 
y + 10y?y2 + 20yíy3 + 15y1y5 + 30y1 ya + 20y2y3 + 24ys 
Qs (Y1, Yo, Ya, Ya) Y5) = - A A A A — ——— A A254A 22222 y 
120 
By Lemma 2.6, we have the relations 
ki=i km-1=1 
Pan (D1)--->Pm) y Thy y Th * NR Thom (2.24) 
ki=1  ko=1 km=1 
n k1 km-1 
Qm (P1;--->Pm) = > hy e Me (2.25) 
ki=1 ko=1 km=1 


Hence, according to the definitions of Stirling numbers of the first kind s(n, k) and Bell numbers 
Y, (n), we can get the well known identities 


s(n,k) =Pr-1 (His a Eh =- iba (9) l (2.26) 
Y, (n) = kIQ4 En Ae HA) = kiC* (101) (2.27) 


which were proved in [49]. 
Therefore, by (2.26), then 


y A A (p+1, (943), (2.28) 


nin? nink* 


The two series are equal, which is an instance of duality for multiple zeta values. (See also the 
discussion on pp. 526-7 of Chen [16]). 
Using Theorem 2.7 and formulas (2.24) and (2.25), we are able to obtain many rela- 
1 
tions involving multiple zeta-star values and zeta values. For example, taking 1, = e 
1 
1 (¿=1,--- ,n) in Theorem 2.7, (2.24) and (2.25), then letting n —> 00, we deduce the results 


(yr m-—1 
Cc Upy”) me 2! De (toy )e C (pm — pi), 
1 m-—1 
iO 2. ds (ty Ya C (pm — pi). 
Setting 1, = a and m = 6 in (2.24), then letting n > oo yields 
20 1 3 1 1 
(109) =- Ec (6) + 3,07 (3) + EC (5) m2 — 56(2)6(3)Im2 + 5,C(4) Ia 
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1 1 
a In?9 — — 2 4:6 
+ 31508) n e las 24 Sa, 


Now we establish some relations between zeta values and series involving Stirling numbers, 
harmonic numbers and Bell numbers. 


Theorem 2.8 ( /49/) For positive integers k and p, we have 


Es s(n, k) Y, (n) k+p 
A =0( a )otro+n, (2.29) 
E S (n, k) (p+1D) — = S (n, p) (k+1) 
2. nin q 2 nin Eno Ea 
a k) Y, 
A Ca 2» ct (nm) _ AE BEDE6 E (2.31) 
n=p n=k 


Proof. Multiplying (2.5) by 27*In”(1 — x2) and integrating over (0,1), using (2.20), we deduce 


that , 
/ InP+* (1 — EN _ (Ap 3 s(n,k) 1 (n) 


z nin 
0 n=k 


By Lemma 2.2, the integral on the left hand side is equal to 


PEA a = (Pp + IC (pH kh + 1). 
0 


Summing these two contributions yields the result (2.29). 
To prove (2.30), we need the well known identity 


1 
Qe recto +1 HD), pen, (2.32) 
0 


Replacing x by t in (2.5), then dividing (2.5) by t and integrating over the interval (0, x), the 
result is 


In*(1.= 
ppal e O a E E wEl (2.33) 
0 


nin 
n=k 


lp? 
Multiplying (2.33) by E Z and integrating over the interval (0, 1), by (2.32), we obtain 
—z 


1 


lava fla (12) O SE 
——=— dltdx = (> pl AA — +0) 
Jj] dida = (1) oa E (<(p+ 1) - HÉ ll (2.34) 


0 0 


Noting that the integral on the left hand side of this equation can be written as 


1 1 1 
P a o P(]— P(i— 
e (1 ne sega 1-u Pp (1 DS => (1 Y ar 
lx t t=l1-y u 1 zx 1 
0 0 


u 0 zx 
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By integration by parts, we see that 


FP (1 In*+ / Pak 
pl 
z 1 1-—¿ 

0 


z 0 í 


dt 


t 


] 
a 
— 
[a] 
3 
AS 
E 
| 
+ 
Y 


1 1 


0 0 0 


=(1+ptk! > S M2) (< (+1) ED) 


Then with the help of formula (2.28) we easily obtain the result. 
To prove identity (2.31), we consider the following integral 


zx 


la co , 
de (1 — x) 10? (1 — (121 may — = (—1Pp! y S Qu) Jo (1 — x)dx 
A 0 


eE rt y RR (n, p) Y, (n 


nn 


Xx 


l zx 
P(i- E k (1 — me 
a E 0 2” pa (1 2" DE (1 pa (1 o tdo 
0 0 0 


= (FP Rtpic (k +1) E (p+ 1) caera A A 


We thus arrive at (2.31) to complete the proof of Theorem 2.8. 


Theorem 2.9 For positive integer m, then 


Y) A — (2, (11) 102€ (2, LUJon-1) 112 (m+ 1)6(8, (1), 


nn! 
n=m 


Proof. Using integration by parts we deduce 


n 


0 7 


Taking x =—1 in (2.36) yields 
1 ya 
Ja In (1 + xjdx = (1) [(-D)”"1n2 — In2+4+ H,). 


n 


0 


Next, we consider the following integral 


1 
InP+1 (1 +) 2 (- 
TA a = (MY 

] on dx =(-1)"m paa n,m 


0 


17 


se 1 q) 
n—-1 n 
In (1 — t)dt = — In (1 — => — — In (1 — =1< 1. 
/ t n( t)d a all %) : n( As x< 


l 
a 
Rl|- 
Mi 
<| + 
a 
X 
S— 

=] 
3 
pp 
Si 
| 
== 
a 
—h 
+ 
A, 
ll 
115 
y |8 
ÁS 
=] 
E 
a 
Xy 
a 


(2.35) 


(2.36) 


= (-1)"m! ln ¿y 19 (mim2 Y cm) 
mE Ec | 


By Lemma 2.2 and Lemma 2.3 


1 00 
pe (1 +2), = ("mi y Am q = (-1)P mi (O, LDm-1), 
0 


con = (19% Y 227 yr (2, (LJ m1), 


Thus, combining the above identities, we complete the proof of Theorem 2.9. 
By Lemma 2.3 and Lemma 2.4, we can get the well known alternating multiple zeta values ( [9]) 


pa” 


SOI E 
a (a) 


Theorem 2.10 For positive integer m, then 


y EA (A ml (3) (2.37) 
n=1 

2 Ym- (n) HE, n=1 __ 1 1 

y a L=m (s (m +1) — Lim+1 (3) , (2.38) 


Proof. To prove the first and second identities. Using (2.20), by a direct calculation, it is easy 
to see that 


il Só 1 00 
PEI aya fer je (ay y a) 
a 0 n=1 


n=1 


Applying the change of variable t > 1 — x to integral above and using (2.11), we obtain 


ln” (1-2), 1 
= (-1"mtbima (2). 
ESTAS Y" milan (5) 


0 


Thus, we deduce (2.37). Using the Cauchy product formula of power series, we can find that 


In(1-x) Y a 
A A SN 2.41 
l+x 2 "Ena", ES (2.41) 
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Therefore, taking into account formulas (2.20) and (2.41), we obtain 


n=1 
= (yr A y 
n=1 


Then with the help of formulas (2.13) and (2.40) we may deduce the result (2.38). 
To prove the third identity. By using (2.5), we have 


1 
n 


1 
mtil=2)... m1 2. s(n,m) zx 
2 dx = (-1) mi Y, e Fe 
0 li 0 


[0,0] Ha 
= (ym! Y A + In "+9, 


Comparing (2.40) with (2.42) completes the proof of (2.39). 


2.4 Some identities on Euler sums 


2. Hn aci C()= 149 
y (=1) RÁ 


e ea a 7 1 j 
SEAT (yal = 26 (8) — 26 (2) Im2+ In82, 
n 8 2 6 


n=1 
2 H,Hn 43 3 1 1 
2 — —=((4) + Éc(2)In22 — =1m%2 — 3Lia | = 
2, 2 A sta (3) 
y H3 43H, H0 + 2H y Li (3) 
——  _— ————J —————————— 0 0 a E 
e 2 
n=1 
00 53 174 AQ 94,7 l 2 
ETA PAR yo nto 2 (31240 (0)+ e (9) 1622, 
n=1 se ñ ; j : 
se 2 2 y 
n M y a 2 
n=1 
SÍ 2-40 Ed SS e So. s 
A n n aL E a Hn n— 
pl a PE e E 
n=1 n=1 . E 
= 2Lia (3) + ina, 
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(2.42) 


(2.43) 


(2.44) 


(2.45) 


(2.46) 


(2.47) 


(2.48) 


(2.49) 


n 


o (43+38, 80 +28) A, 
Za (ei ll ME = 24 (< (5) — Lis (3) ¡ (2.50) 


2 
o Hi4+ 84,40 + 68280 +3(447) +6H0 
A aa a = 24Li, (5). (2.51) 
o (42-HP) A, eN 
E > 
1 
= 12 (4)In2- 3 324 —In%24+ = S ia [|< 
C (4) In mios qn 2+ (3) in 2 + 18Lis (5)m 
+ 24Li, (3) —24c (5), (2.52) 
a a _ 2 
o Hi+8 4, HO + GEO +3 (E) +6H0 
A (aL) 
= 56 (5) +30 (2) € (3) + FC (9) 124 3€ (3) In%2 + 2€ (2) ln%2 + ¿1m02, (2.53) 
00 (413 — 3H. O 0) A, E 00 (4 — 19) ña a 
y A Y) 22) 
n=1 n=1 
a HnHn m—1 2 Ha n—1 : 1 Ll 5 
AA o O o A : 
jos (1) 2 ) óLñs (5) a (2.54) 


These identities can be obtained from the main theorems which are presented in the section 2.3. 

For example, setting m = 1,2,3,4 in (2.37) respectively yields (2.43), (2.44), (2.46) and (2.51). 

AA 
k 


the results (2.47) and (2.53). etc 


Taking 2. = and m = 4,5 in (2.25) and Theorem 2.7, then letting n > 00, we obtain 


3 Some explicit evaluation of Euler sums 

In the section, by constructing the integrals of polylogarithm functions, we establish some rela- 
tions involving two or more Euler sums of the same weight. Then we use the relations obtained 
to evaluate alternating Euler sums, and give explicit formulas. First, we need the following 


lemma. 


Lemma 3.1 ( /48/) For positive integers n, m and real x € [—-1,1), we have 


n—1 = he (—D”" (Ex 
Join 0d Y y imei) E EF 
0 ¿=1 k=1 
cp A (1-2) (a” -1). (3.1) 
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3.1 Main theorems and corollaries 


Theorem 3.2 Let m,p be integers with m,p > 0, |x|, lyl,|z] € (0,1), we have 


> Li (62) Gn (05 4) — Lip (Y) En (M3 2) 


n 
n=1 
= Li, (y) 2 ma Li (06) 2, a 
+ Lim (2) Lip+1 (24) — Lip (y) Lim+1 (22), (3.2) 


where the partial sum C, (m; a) is defined by 


n A 
ar 


En (m; a) = m7 (a € C). 
k=1 


By convenience, we let (o (m;a) =0. It is obvious that 


Calm l)= am and Cn(m;-1) = BN. 


Proof. Let us consider the N-th partial sum 


No. 
Sy (2,9, 2) = Y Lim (4) 6n (Pi) — Hip (9) en (132) 
n=1 
il N 
= Lim (2) y, ón 2 y) Z" — Li, (y) a ón un z) n 
ml n=1 
N n yen N n gh gn 
= Lim (x) y e nkP — Li, (y) y y, nkm 
n=1 k=1 n=1b=i 
N N yz? N N kn 
= Lim (2) 7) — Lip (9) 70 
k=1n=k n=T—k 
sd N 
= Lim (2) ) vilo) ia, — Lip (y) Y evi) eri) 
k=1 a 
= Un (1; 2) (Lim (12) Cv (py) — Lip (y) Cn (m; 2) 
N N 
+ Lip ( e al 2) a Y al 2) 
4 h=1 


Then letting N > oo in above equation, we get the required identity (3.1). The proof of Theorem 
3.2 is thus completed. 
Taking (x,y,2) > (1,1,1), 1,—1,1) and (-1,—1,1) in Theorem 3.2, we can get the following 
corollary. 


Corollary 3.3 For positive integers m and p, then 


00 (p) (m) 
> E? AS" 
n=1 


n 
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n=1 de 
+5(m)¿(p+1)-C()c(m+1), (m>2) 
2 E LOAPD Ema - l 
A 
n=1 


+€(m)E(p+1) —C(p)C(m+1). 


The main Theorems of this section can be stated as follows. 


Theorem 3.4 For positive integers m,p and real a,b, x with the conditions ja] = 


lz] < Jal? or [bJ7*, we have 
Gn (1; az) Cn (mi b) moy" 6n (1; b2) Gn (p; a) 
Y arar? El 1) e pnpm 


1 


p-1 00 . m= bs _ 
= 2D Liga (am) Y, en ma) - Y (1) Lim+1-: (b2) Y - le : 


¿=1 n=1 i=1 n=1 
2 n :D n —n ya = 

+ (Ya) YD -a 7) -(-)"In( 100) 0 
n=1 n= 


(3.5) 


lb] > 1 and 


qe 


(ar =p"). 


(3.6) 


Proof. Applying the definition of polylogarithm function Li,(x) and using the Cauchy product 


of power series, we can find that 


mula 


- Yanira) tela 


Then using above identity it is easily see that 


z 


Lip (at) bel bt) m=lr: qm = E n—1x ; 
j =D Ed ye Gn(p; a rf Lisa (dEJdt= 2 En (mi; b) q Li, (at)dt. 
0 me 


o 


0 
Then with the help of formula (3.1) and the elementary transform 


z bx 


y SL d= a 1 de 


0 0 


by a direct calculation, we may deduce the result. 


Theorem 3.5 For positive integers m,p and real a,b, x with the conditions la] = 


lx] < Jal"? or [bJ7*, then the following identity holds 


e A Pa 


ar? nm 
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lb] > 1 and 


p-1 p-i m-—1m-i 


= Y (- D"Tipy2- 5 (ax) Lig 45m (5, b) =Z 20 VA TLim+9-¿ 5 (bx) Li;y ;p (2, a) 


1=1 ¿=1 =1.y=1 
+ (-1)p1n (1 — az) (Lip41,m (2,b) — Lip41m e Ze )) 
— (1) mln (1 — bx) (Lim+1p (2, a) — Lim pia) 


Gn (1; a2) ón (m;b) | hb) (pa) 
ty anyo+l "my jnaqm+l 
+ (-1 Lia (02) Lipm (a7*,b) — (-1)” Liz do Dip (00), (3.7) 


where Lis « (1, y) is defined by the double sum 
Lx 
Line (00) 25 [| < 1,1u] >0,|ay| <1,5,t >0 


Of course, if s > 1, then we can allow |xy| = 1. 


Proof. Replacing x by t in (3.6), then dividing it by t and integrating over the interval (0, 2) 
with the help of formulae (3.1) and 


] m0 = —Lio (ax), 
0 


by a simple calculation, we have the result. 
Letting (a,b, x) > (-1,1,1),(-1,—1,1) and (p,a,b,x) > (1, —1,1,—1), (1, —1, —1,—1) in The- 
orem 3.4, then with the help of Corollary 3.3 we obtain the following corollaries. 


Corollary 3.6 For positive integers p > 2 and m > 2, then 


m-—1 p-1 
CD" Si5m + (D"Smip = Y (DC (m+1-i) ple p+1-4)Sm, 
1=2 i= 
+€ (9) Sim — E (m) Sip +5 (m)5 (p+1) -C(p)C(m+1), 
+ (-1Y (Sm,p + Sm,p)1n 2, (3.8) 
m-—1 p-1 
(—1PSimp => (—D)" Sig m = (Dc (m ap ¿) Sp SN (=D (p plo ¿) Smi 


2 1=2 
+6 (9) Sim — € (m) 815 +€(m) € (p+1) — E (p)C(m+1) 
+ (-1P (Smp + Sm,p) In 2 — (= ÓN (Sam + Som) ln 2. (3.9) 


Corollary 3.7 For integer m > 1, then 


m-—1 

(DMA Sia mm + Sima = Y) ED E (m4 1) 877 + (1977 (Sim + Sim) In2, (3.10) 
i=1 

CDO"S im + Simi => EDO C(m+1i)S5,. (3.11) 
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Similarly, letting (a,b,x) > (1,1,1) and (-1,—1,—1) in Theorem 3.5, we can get the 


following two corollaries. 


Corollary 3.8 For positive integers m > 1 and p > 1, we have 


p-1 p-1 
(-1P Som,p — (=1 Sm =D) (- El (p+2—i— jSmitj 
i=1 j=1 
m=1m-—i 
Y Pe m+ 2 id) 
i=1 j=1 


- (-1PPSimp+1 + (-D)"mSipm+1 
+ (DC (2) Smp — (1302) Spym- 
Corollary 3.9 For positive integers m and p, then 


m=—1m-—1 
(—1D)" Sopm 1*S2mp.= = y y ¡ica (m FL=4 9) Sam 


i=1 j=1 
II A 
pos 


i=1 j 
e CO" msm. = COS msn 
APT SEAS 


Theorem 3.10 For positive integer k and 1< í1 < i2 < -*- < iy with x; € [-1, 
following relation holds: 


k k-1 3% 
E e 
pr JAI : TR 


j=0 1<h <-<i<k y += (1, ++4,) 


where by convention, the term corresponding to y =0 is A (1 + +idp 10007 


Proof. Consider the generating function 


(3.12) 

(3.13) 
1], then the 
E (3.14) 
ú, 


Fr Es BE a 1 dy; 5) — Cn (1 +: +dpj 0100 * Tp) E, (x € (-1,1)), 


n=1 
where Xy = (21,-:* ,Tp). 


Using the definition of the numbers C,, (m; a), we have 


k 
(l=x)Fi(u Xp) = nl) = bn i++] [ 2; qn 


k k k 

E E 
11 E 1005) = Chcd ) 1d | | Ds zx 
j=1 j=l j=1 
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k 
II 5 


00 k k n Fl 
=> lots) - (tre) -E) - E a 


E 


Us 


Tota) Te 


[0.0] 
DE OS + ¿Y Y. RR a" 
1 


n= 


Multiplying this equation (3.15) 
identifying the coefficients of x” 


, we may deduce the desired result. 


j=0 1<h<--<l¿<k y cea (6 ++) 


(3.15) 


by 1/(1 — 2), using the Cauchy product of power series, and 


3.2 Identities for alternating Euler sums 


It is clear that from Theorem 3.4, 3.5, 3.10 and Corollary 3.6-3.9, we can establish many identities 
involving two or more Euler sums of the same weight. Some examples on alternating Euler sums 


with weight < 5 follows: 


= n ¿Ml n— = Hr qe n— 
NN Pray y A yt Be (4) + 3600) 102, (316) 
n=1 n=1 
> H, HO E? 7 
Y = Ey - y ON ES (4) +68) n2— =e (2) In?2, (3.17) 
n=1 n=1 
2 (27 if SHA, 43 15 35 

n—1 n n__ 10 99 
2 A e pS A GOO) CS) + 0 (4) In2, (3.18) 

n-1 nin n=1 e: 9 
2 n A 22, SED = TAS! ¿(20 (3) + 754) In2, (3.19) 
2. HH o E Ñ 
ON Haya Y Aaa 0 (3) 1224 202) C (8) — FE (4) m2, (320) 
n=1 n=1 
00 (2) 00 25 
y aL +) E = ALia (3) In2+ ina dh Lot ) In22 — € (2) In32 
n=1 n=1 

5 23 
«E o ES ÓN In 2, (3.21) 


3 49 
Al 16* (2c(3) + 6 (4) In 2, (3.22) 
= —2Lia (3) In 2 — In + —£ (3) In?2 + <C (2) In*2 
+ 5502) 0(3) — (9) In2, (3.23) 


(2 5 o (2) 7) 00 7 (2) 
Hr Hn pe Hñ Hp PE Han Hn ez 
y A UY PA E 
n=1 n=1 


= —4Lis (3) — yin + EG (5) + le (3)1*2= 5 (2) In%2 + En (2) 6 (3) — 3€ (4) In 2. 
(3.24) 


These identities can be obtained from the Theorem 3.4, 3.5, 3.10 and Corollary 3.6-3.9 which are 
presented in the section 3.1. For example, taking k = 3, (i,,12,13) = (2,2,1) and (21,12, 13) = 
(1,1, —1) or (—1,—1,—1) in Theorem 3.10, then letting n > oo yields the formulas (3.18) and 
(3.19). 

In (2.5), setting k = 3,4 and x= —1, using (2.18) we can get the following two identities 


H3 3H, HO + 247 1 21 1 
y ae Ml o MN A = 60 (4) + ze (2 1n%9 = 31092 - 75 (3) In2— 6Lia (3) 


n 

n=1 
(3.25) 

2 
o H14+8H,HÓ — 68280 + 3 (44) 6H 
A o A 
n 
n=1 
e 1 e ; ; 
= 24Li; > + 24Lia > In 2 + In"2 + e (3) In*2 — 24c (5) — 4€ (2) In*2. (3.26) 


3.3 Closed form of alternating Euler sums 


In this subsection, we use our equation system to obtain explicit evaluations for alternating 
Euler sums with weight < 5. 

Note that the linear sums Sy, Sp, and S5,¿ are reducible to zeta values and polylogarithms 
when p+q is odd (see [12,20,23,45]). The quadratic sums Si2 357, S12375), 9113759 911,2m> 912 2m 
can be computed by the paper [48]. The quadratic and cubic sums Sj2 3, S72 3, S12 3, S13 3, Sisa 
and Sis ¿ can be evaluated by Bailey, Borwein and Girgensohn [3], we get 


[0,0] 
He í 1 1 21 

y 1) = 6Lis (5) + 6Lia (5) 1m2+=HIn92 + << (3) In?2 
n2 2 2 5 8 


n=1 
9 e 

- 355) - €(2) In%2 — Ec (2) 5 (3), (3.27) 
E H? n-1 Ñ 1 A 1 2 5 7 2 
2 ya EU) = 4Lis (3) + 4Liz (3) In 2 + 75 2+ ri (3) In“2 

55) 55 01002 ECC), (3.28) 
RR .f1 Lua E 2. 167 1 3 
Dl — ALis (3) e 2 +35 (3) In 2-37 (8)+ 3002) In > 

+30) C(8) + EC (4)1m3, (3.29) 


2 a (Ms 
y —(-1)"* = —4ALi, (3) In 2 — qa 2-— 3266) + ( (2) In*2 
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3 19 
+2c(2)c (8) + Ec (4)1n2 (3.30) 
Hence, from identities (2.46)-(2.49), (3.16), (3.17) and (3.25), we deduce the following six closed 


form 


00 25 

Hr Hn n—1 43 1 2 1 4 1 1 
A a E q Lo E a A a E ON 31 
E E E E aa (3). si 
o Ada n—1 93 1 4 1 2 ] ! 
E E (3) di 
2 A e el 59 13 e. A al 
Ll (3) A AE ens 
y HH 12 = oia (2) 4 E1mt24%e (3) m2 — Le (2) 1022 (4) (3.35) 
ES > TS E 
00 3 
$ ma (pde = (4) + a (2) In?2 — 212 ES ze (3) In2. (3.36) 
n=1 


Note that formula (3.36) was also proved in [33] by another method. 
In (2.5) and (3.8) of [48], letting s = 3 and (p,!) = (1,3), respectively, and combining (3.28)- 
(3.30) of this paper, we give 


2. HaHn Pr A 7 5 
+ =C (2) In%2 + 4€ (4) In 2 — ze (2 (3.37) 
WA Hana o. (1 1.5, 193 1 2 
2 —— =2Lis (3) = 02 — E (5) — ¿E (8) In?2 
+ e (2) Inó2 + 4 (4) In2 + SE (2 (3). (3.38) 


Combining (3.20) with (3.37), and (3.23) with (3.29), we have the two results 


A >. 37 A 
2 n [AJA =2Dk (3) - e a 0 (3) In%2 — 150) JE 560) 1n39 
+6 (2) € (8) + E (4) In 2, (3.39) 
a (3) Ls (5) ln2— zp In 24 30 (3) m2 + 37 (5) 
1 3 1 49 
+ 562) 1092 — 502) (3) — ¿(9 In2. (3.40) 
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By (2.18), letting m = 4, z = 1, we arrive at the conclusion that 


1 
In? E a H3-3H HP +2 EH 1 
1 AC 
(n +1) 
E (1 ds 21 a , 
= —24Lis a 24Lia a ln 2 — ¿ 2— E (3) 1n%2 + 24€ (5) + 4€ (2) In?2. 
(3.41) 
Substituting (3.27) into (3.41) yields 
3 HH 1)? = -4Li aaa (o da a Te(3)1 ea 
La AO E A 
23 2 15 
+ =C (5) + 50 (2) 1n%2 + =C (2) € (3). (3.42) 
8 3 16 
Combining (3.22) with (3.42) gives 
SE ay, a =<(3) Ia 23 (5) + Lc (9) In32 
Lona PAX) 30 ARE go 
3 49 
— 75 (2) (3) + E (4) In 2. (3.43) 
Theorem 3.11 For positive integer m, then the integrals 
Fln2Ma in? (1 Fln2Madn? (1 
PP cla” (1-2) ,. na [2 aa (1+2),. 
l+zx lx 
0 
are expressible in terms ofln2, zeta values and linear sums. 
Proof. Using (2.32) and the well known identity 
Aloe NN g(m+1) 
] a dx = (-1) mi (E(m+1) — HN a men, 
0 
we can find that 
1 
ln” gin? (1 — 2) > HEN 2 A a 
o A AO Y AN AS A les - 
] o 2A-)"m1 4) a 5 Po 1) 
0 n=1 m=l 
+ (-1)”m!€ (m +1) In?2, (3.44) 
1 
ln”xIn? (1 + 2) de 2 H, EmeD E green aia 
0 n=1 n=1 
+ (1) mic (m + 1) In?2. (3.45) 


From Corollary 3.7, we see that the quadratic sums 51(2) 1 and 9,5, 1 are reducible to linear 
sums and zeta values. Thus, Theorem holds. 
Setting m = 2 in (3.44) and (3.45), and combining (3.39) and (3.40), we obtain 


1 
PEA $e (3) _ nó +00) + 5 (2) m2 — << (2) 6 (3) +C (4) m2, 


l+zx 2 
(3.46) 
1 
In?zIn? (1 + 2) dl ol A 9 47 
2 13 49 
+ AS (2) Iné2 — eN (2) € (3) — e (4) In 2. (3.47) 
Now, we consider the generating function 
= E, (-1<zx<l). 
n=1 
By a direct calculation, it is easy to see that 
flo (1—t) A yr 
=2 2 |) =———dt =>. 3.48 
y 0y+2] me pe (3.48) 
0 ut 
By integrating, we may rewrite (3.48) as 
2 qe 2n(1-+t 2102 (1-4 
y En Ba nt 3 (1 ») | A 2] ESO (3.49) 
A=1 n+ 1 0 1 + t 0 1 + t 
1 
Multiplying (3.49) by == end integrating over the interval (0,1). The result is 
y 
2 HA 5 15 9 12 2 
E E A a 0) E) pa) Me A 
A + Ma Al E? Dr 
1 
= Hno Hon 1 f In?xlin? (1 — 
pe y - A ] fala? (12). (3.50) 
n 2 l+zx 
n=1 0 
Substituting (3.29), (3.37) and (3.46) into (3.50), we have 
NE 
HO E a E e UN 7 29 
2 (17 = —4Lis | 5) + 1092 + -C (3) In%2 - =C (5 
2 io (5) + qna (0% elo 
1 15 23 
- q (2) Inó2 + ra (2) € (3) — ra (4) In 2. (30510) 


Inz 
Applying the same argument as above, replacing x by —= in (3.49), then dividing this by al 
z 


and integrating over the interval (0, 1), we arrive at the conclusion that 


o $ 4) 
Alpaca = (mo 2 
2 aLis (5) + 105 Pg 
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Letting (a, ,b,p,m) = (1,1,-1,2,1), (-1,-1,1,1,2), (-1,1,—1,1,2) and (-1,-1,-1, 


y 5 (2) In32 + == DUS + Ec (4) In2. (3.52) 


a 
Theorem 3.5, and combining identities (2.51), (3.18), (3.19), (3.21), (3.26), (3.37), (3.38), (3.40), 


(3.42), (3.43), (3.51) and (3.52), we obtain the following eight evaluations 


00 py (2) 
Y) TL 0) ESC), (3.53) 
m=1 
o A 3 25 
2 DA 66) + 352) 0 (3) — 755 (4) In2, (3.54) 
2. HH? o cl 2 29 7 
Za a ALis (3) + 4Lia (3) m2 502 + o (5) + ni (3) In?2 
0% 0 (2)C (8), (3.55) 
00 (2) 
a (1)? = —2Lia (3) In2— In = E (3) In*2 + 5 (2) In%2 
m=1 
7 
+30) 0(3) + 0(9) 102, (3.56) 


25 
> An An 0 Sl E oa 2 1 ñ 


+55 (8) - Lc (4) m2, (3.57) 
o (4 y 
o E = 16Lis (3) + 8Lia (3) In 2 + Ina - € (5) — 50 (2) In%2 
n=1 
+ 20(2)6(3) + HC (4) 1m2, (3.58) 


o (qoy 1 1 1 299 2 
y Eo = —16Lis (5) = 8Lia (3) ln 2 — 51092 + (5) + ES (2) In%2 
Í 39 
=350)5(8)= 75 (4) m2, (3.59) 
o rá 
y Anaya = 8Lis (3) + 4Lia (3) In2 + Im + Ze (3) In?2 — 75 (5) 


Ñ e (2) Ins9 — Et (203) + =< (4) In2. (3.60) 


In many other cases we are also able to obtain a formula for other Euler sum of weight five 
explicitly in terms of values of zeta, logarithm and polylogarithm functions. For example, we 
also use the method of the present paper to obtain explicit evaluations for the following sums 


9) (S 5 
Hr Hna mi (1 Los 3 9 193 1 3 
————(-1 = —2L —- —1n*2 — — 1n%2 — — — —( (2) In?2 
2 E is (3) + 2 ¿E (8) In (50) ¿6 0) In 


30 


11 5 
+ (9 In2- ¿6(2)0 (8), (3.61) 
= HH mi (1 115 3 9 19 1 3 
2, mn (=1 = 2Li, (3) In 2 + 371 2+ 36 (3) In*2 — 33% (5) — a (2) In?2 
11 1 
+5 (9In2+76(2)0(8), (3.62) 
SH o 1 (1 1 165 1 
2 nz = 4Lis (3) + 2Lia (3) ln2+ 32 - EE (5) — ES (2) In?2 
+00) 0(3) + 775 (9) In2, (3.63 
ZHHn, mnmt o: (1 E 2. 197 1 y 
2 7 (-1) = 2Lis (5) — q 2 — 36 (8) In%2 — an: (5) + ES (2) In?2 
3 109 
32008) + 37564) In2, (3.64) 
BB e 2. 899 2 , 
2 (DD =14Lis (3) — ¿q e2— ¿E (8) In*2 — EC (5) — E (2) In?2 
1 193 
+7520)08) +35 m2, (3.65) 
00 pr(2) 57 
Hi (== 6 (3) In2— 1n59 — € (2) In%2 + a (2) € (3) 
a. E 2 4 8 
+ e (4) In 2, (3.66) 
AA Dis 9 A 7 A 
2 (== o (3) — ea — 756 (3) 1092 + 750 (5) + 3762) In*2 
1 37 
= 5 (M0(8)+ GC (4) In2, (3.67) 
SH HN E 1 E zo 155 1 , 
2 ==, ED =-6Lis (5) + q 10924 50 (3) 1924 33€ (5) — 77€ (2) m2 
+ Lc (2)c(8) — 2 (4) m2. (3.68) 


Some of the above identities can be found in [45]. Moreover, we use Mathematica tool to check 
numerically each of the specific identities listed. The numeric results prove that the obtained 
formulas are correct. 


3.4 Some evaluation of Euler-type sums 


In this subsection, we will establish some explicit relationships which involve Euler type sums 
of the form 


Smima..mp(1 521) +00,2p,2) = Y de p,k,r EN 
n=1 


(n +7) (n + k) 
and (alternating) harmonic numbers, where x, a; = +1 and r < k. We develops an approach to 
evaluation of Euler-type sums Sm,ma...mp (Pibes ¿E x). The approach is based on integral 
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computations. 
Theorem 3.12 For integers 0 <r < k and i1 < ía < -**: < ip, then the following relation holds: 


(kr) (Sas (r,k;x1,-** ,tp;1) — O miidta] o 1)) 


Ñ > 5 apa 
y ¡=1 2, at (a+ +11) 
J 
kr | p-1 IT Cn (i1,; 21.) II za 
a YN Y Y — E (3.69) 
i=1 | ¡=01<t1<--<!¡<p y e+Hip= (6, +=+u,) (n+r+i-1) II xr 
1 
b=1 ? 


Proof. From (3.15), we have 


¡0.0 


Pp 
[O (2) — (1 + + ip201 00 2p) pu" 


1% eta al 

oo 1 Jj p-1 Il Ón Uan.) ll. La, qa 

== IZ] p-j r= a= 

e = y y 1+ Hip + > ) , (=1) == A == . eS (3.70) 


n=1 j=0 1<11<--<lj<p a O) 


Multiplying (3.70) by x"7! — x*=1 and integrating over (0,1), by a simple calculation we may 
easily deduce the desired result. 
Hence, from Theorem 3.12, we can get the following corollaries. 


Corollary 3.13 For positive integers m1, ma, k and r > 0, we have 


00 gm gl) gira 


(=") 2, (n ++") (n+k) 


n=1 
Ñ k—r > go A am Ñ 5 e 
=== neln+r+j-1) nu(n+r+j-1) numa (n+r+j-1) |” ] 

oo. (mi) (ma) (ma +ma) 

n Hr = Há 

(lo 0) Y 

n=1 (n Al r) (n + k) 

k—r e 


00 arta aria 9 
a a ii 
neln+r+j-="1) ne(in+r+j-1) nutme(n+r+j—1) 


gm) (ma) great 


n n 


(n ++") (n+k) 


k=r 00 arica al Ns 
2 e E Su E 
, naln+r+j3-=1) nu/in+r+j—1) nutre (n+r+j—1) 
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where the Euler-type sums > d oy PE (p,m,j € N) can be evaluated by 


rn? (n E 3) 
the following formula 


Note that the closed forms of the following Euler-type sums can be found in [39, 53, 54]. 
9) gl o) (m) 
Hr 
—_—_——————— k. 
2 (n+k)” EE (n+"r) (n+k) a 


Corollary 3.14 For integers m; > 0,ma > 0,m3 > 0 and k > r > 0, we have 


ES aga gin) plmi+ma+ms) 


da (n+r)(n+k) 
kr oo AMD y (m2) gim ln) gl gln) 
A 
E Bl Bl» aga 


It is possible that some of other Euler type sums can be obtained using techniques of the 


present paper. 
Acknowledgments. The authors would like to thank the anonymous referee for his/her 
helpful comments, which improve the presentation of the paper. 
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